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Abstract. We consider an initial- and Dirichlot boundary- value problem for a fourth-order linear 

1 I , stochastic parabolic equation, in two or three space dimensions, forced by an additive space-time white 
1 noise. Discretizing the space-time white noise a modeling error is introduced and a regularized fourth- 
• order linear stochastic parabolic problem is obtained. Fully-discrete approximations to the solution 

( ' of the regularized problem are constructed by using, for discretization in space, a standard Galerkin 

(-H ^ finite element method based on piecewise polynomials, and, for time-stepping, the Backward Euler 

. method. We derive strong a priori estimates for the modeling error and for the approximation error to 

the solution of the regularized problem. 

1. Introduction 

T— I ■ 

> ■ 1.1. The main problem. Let d = 2 or 3, T > 0, = (0,1)'' C and {n,J^,P) be a complete 

OO , probability space. Then wc consider an initial- and Dirichlet boundary- value problem for a fourth- 

(N . order linear stochastic parabolic equation formulated, typically, as follows: find a stochastic function 

^ ; u:[0,T]xD^R such that 

\0 ' dtu + A^u^W{t,x) \/{t,x) e {0,T]x D, 

(1-1) A'"^.(t,.) 1,^-0 V<G(0,T], m = 0,l, 

O' u(0,x) = \/xeD, 

• T^ , a.s. in il, where W denotes a space-time white noise on [0,r] x D (see, e.g., [H]). The mild solution 

^ ' of the problem above (cf. [S], [T^), known as 'stochastic convolution', is given by 

■ — ' (1.2) u{t,x)= f f G{t~s-x,y)dW{s,y). 

Jo J D 

Here, G{t; x, y) is the space-time Green kernel of the corresponding deterministic parabolic problem: find 
a deterministic function w : [0, T] x £) ^ M such that 

^tW-fA^tw^O V(t,a;) e (0,r] X D, 

(1.3) A™u;(t,-)|,„ =0 Vfe(0,T], m = 0, 1, 

w{0,x) — wo{x) y X e D, 

where wq is a deterministic initial condition. In particular, we have 

(1.4) w{t,x) ^ [ Git;x,y)wo{y)dy \f it,x) £ iO,T] xD 
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and 

(1.5) G{t; x,y)^ ^ e"^' * e^{x) e„(2/) V {t, x, y) e{0,T]xDx D, 



where tt^ |a|^d, laln-i ■= {j2i=i Q^f) ^ and £0(2) := 22 HiLi sin(Q;i ttz,) for all z e D and a £ N''. 

1.2. The regularized problem. Following the approach for a second order one dimensional stochastic 
parabolic equation with additive space-time white noise proposed in [I], we construct below an approxi- 
mate initial and boundary value problem: 

For N^,, Ji, G N, define the mesh-lengths At :— Ax := j-, and the nodes t„ nAt for 
n = 0, . . . , Ni, and Xj :— j Ax for j = Q, . . . , J^,. Then, we define the sets AC {1, . . . , A^*}, 



Mi 



{1, . . . , J4, T„ := (t„_i, <„) for n £ AC, Dj := (x^-i, x^) for j G J-*, I?,, Jlti 
for /X e i//, and Sn,fj, ■= Tn x for n G AC and fi £ J'^. Next, consider the fourth-order 
linear stochastic parabolic problem: 

dtu + A'^u^W in (0,r] x D, 
(1-6) A"M(i,.)|,^ -0 Vte(0,r], m = 0,l, 

u(0,a;)=0 MxeD, 

a.e. in 17, where 

t?(<,a:):=^^^ E E -^---m ^) ^"'^ V (i, x) € [0, T] x A 



/ IdW ,\/neU*, V/iSj'i^, 

J s„.„ 



and A's is the index function of 5 C [0, T] x D. 
The solution of the problem ()1.6|) . according to the standard theory for parabolic problems (see, e.g, [19]), 
has the integral representation 

(1.7) u{t,x) ^ [ [ G{t- s;x,y)W{s,y)dsdy V {t,x) e [0,T] xD. 



JOJd 

Remark 1.1. The properties of the stochastic integral (see, e.g., [26]j, yield that i?"'*' ^ A/'(0, Ai(Aa;)'^) 
for all {n, pt) e AC x . Also, we observe that E[i?"'^ i?" ] = for (n, /i) ^ (n', /i'). Thus, the random 
variables (^"''')(n,Ai)eAA* x j<i o'^s independent. 

1.3. Main results of the paper. The problem (jl.ip is a linearized formulation of the stochastic Cahn- 
Hilliard equation (cf. [5], [TU]) which was introduced by Cook [S] for the investigation of phase separation 
in spinodal decomposition (see, e.g., [16j . [12j). For the convergence analysis of approximation methods 
for fourth-order stochastic parabolic problems driven by a space-time white noise, we refer the reader: 
to [4] which considers a finite difference method for the stochastic Cahn-Hilliard equation, and to [22], 
[13] and [14' which consider time-stepping methods for a wide family of evolution problems that includes 
, while the finite element method is not among the space-discretization techniques considered in ^13| 
and [11]. Also, we refer the reader to [5], [I], [17], [H], ^\ and [5] for the analysis of the finite element 
method for second order stochastic parabolic problems. 

In the paper at hand, we extend some of our results for the ID case, given in [18], to the 2D and 
3D case. In particular, we consider approximations of the solution u of (|1.6p produced by the Backward 
Euler time-stepping combined with a C^— finite element method, and analyze its convergence to the mild 
solution of (jl.ip . This error splits in two parts: the modeling error that appears by approximating u 
by u and the numerical approximation error for u. The estimation of the modeling error is achieved, in 
Theorem 13. 1[ by obtaining the inequality 

i_ 

(1.8) Tn&-^[j(j\u{t,x)~u{t,x)\'^dx\dpY<C Ax^'^ + At^ 



Moving to the direction of building approximations of u, we let M £ N, (Tm)m=o ^^"^ nodes of a partition 
of [0,r], i.e. To — 0, Tm = T and r,„_i < r™ for ra — 1,...,M, and define '■— (Tm-i,Tm) and 
:= — Tm-i for m = 1, . . . , M, and set fc^^x := niaxi<m<M km] also, we let M/i C Hq{D) n H^{D) 
be a finite element space consisting of functions which are piecewise polynomials over a partition of 
D in triangles or rectangulars with maximum diameter /i, and define a discrete biharmonic operator 
Bh : Mh ^ M,, by 

(1-9) {BhV,x)o.n^{Av,Ax)o.o "i^.x^Mh, 

and the usual L-^(Z?)— projection operator Ph : L'^{D) — > Af;i by 

{Phf.x)..o = {f.x)..o yx^Mh, yfeL^iD). 

To construct approximations to u, we employ the Backward Euler finite element method which begins 
by setting 

(1.10) 0, 
and, then for m = 1, . . . , M, finds C/™ G Mh such that 



(1.11) c/r - u;^-' + km BhUj:' = / PhW ds. 

Estimating the numerical approximation error for u, we derive first, in Theorem l6.21 the discrete in time 
(Lp(L^)) error estimate: 

(1.12) ||^fc„^^(^y" |C/^(x)-S(T™,x)|'dx^dp| < C [(fc„_)^ +e-U''- 

where: v — v{r, d) is given in (|5.13p and depends on the space dimension d and a parameter r £ {2, 3, 4} 
which is related to the approximation properties of the finite element spaces Mh (see ()2.19p ). Only for 
the needs of the proof of (|1.12p . we introduce a space-discrete approximation Uh of u and analyze its 
convergence in the L'^{Lj,{Ll j) norm (see Theorem 15. 2p . Also, assuming that the nodes {Tm)m=o 
equidistributed with At — km for to = 1, . . . , M, we arrive at the discrete in time L^(Lp(L^)) error 
estimate (see Theorem 16. 5^ : 

(1.13) ^max (^y" \UJ^{x)~u{Tm,x)\^ dx^dPy < C e^^ At^-'' + h"-'^ . 

To get the estimate above, first we consider the Backward-Euler time-discrete approximations of u and 
analyze their convergence in the discrete in time L^(Lp(L^)) norm above (see Theorem 14. 2[) . and then, 
we compare the Backward-Euler fully-discrete with the Backward-Euler time-discrete approximations of 
u (see Proposition 16. 4p . This procedure gives the possibility to estimate separately the space and the 
time discretization error in constrast to the technique used in and [5] for second order problems. The 
reason of including the error estimate l|1.12p in addition to the stronger norm error estimate (|1.13p . is 
that the order of convergence is slightly better and it allows a nonuniform partition of the time interval 
which is non standard among references (cf. [27], [13], [14], [25], [2]). 

We close the section by an overview of the paper. Section [5] introduces notation, and recalls or prove 
several results often used in the paper. Section [3] is dedicated to the estimation of the modeling error. 
Section [5] defines the Backward Euler time-discrete approximations of u and analyzes its convergence. 
Section [5] defines a finite element space-discrete approximation of w and estimates its approximation error. 
Section [6] contains the error analysis for the Backward Euler fully-discrete approximations of u. 

2. Notation and Preliminaries 

2.1. Function spaces and operators. We denote by L^{D) the space of the Lebesgue measurable 
functions which are square integrable on D with respect to Lebesgue's measure dx, provided with the 
standard norm ||(?||o.£3 {J^\g{x)\'^ dx}^ for g G L^{D). The standard inner product in L^{D) that 
produces the norm || • ||o,d is written as (•,-)o,£>, i.e., {91,92)0,0 ■— 51 (x) 52 (cc) dec for gi, 92 £ L'^{D). 



For s G No, H^{D) will be the Sobolev space of functions having generalized derivatives up to order s in 
the space L^{D), and by \\-\\,,n its usual norm, i.e. \\g\\,,D ■= {EagNg, |a|^d<s II^"5I1o,d}' 9 ^ H^iD). 
Also, by Hq{D) we denote the subspace of H^{D) consisting of functions which vanish at the boundary 
dD of D in the sense of trace. We note that in Hq{D) the, well-known, Poincare-Friedrich inequality 
holds, i.e., 

(2.1) II.9II0.Z, < C^^ WVgWo.n yg e H^{D), 

where ||V,;||o,„ := (EaeNg>|,.=i ll5>llo,z,) ' for v G H^iD). 

The sequence of pairs {{Xa,£k)} qgn*^ is a solution to the eigenvalue/eigenfunction problem: find 
nonzero (p G H^{D) n Hq{D) and cr G K such that —Aip = aip in D. Since {Sa)aeN'^ is a complete 
(•, orthonormal system in L^{D), for s e M, a subspace H (D) of L^{D) (see [53]) is defined by 

and provided with the norm ||v||hs := ( X^aeN'^'^a ("^j ^a)o,u ) ^ S H*(L'). Let m G Nq. It is well- 

known (see ^23j) that 

(2.2) H™(L') = {w e ff"(I?) : A*f|o^=0 if < i < f } 
and there exist constants Cm, a and Cm^s such that 

(2.3) an,A \H,„,n < Iklle™ < C^.s \\v\\„..n V« G H"(i?). 

Also, we define on L^{D) the negative norm || • 1|_„,,d by 

:=sup{^^^f^ : ^ G H™(i?) and ^ ^ o} \fveL^{D), 

for which, using (|2.3p . it is easy to conclude that there exists a constant C_„i > such that 

(2.4) IMU„,o < C_™||«||h_„ V«gL2(I?). 

Let L2 = {L'^{D), (•, •)o,d) and £(L2) be the space of linear, bounded operators from L2 to L2. We say 

that, an operator F G £(^2) is Hilbert-Schmidt, when ||F||hs := {X^fe^i llr^fcHo.o}^ < +00, where ||F||hs 
is the so called Hilbert-Schmidt norm of F. We note that the quantity ||F||hs does not change when we 
replace {efcj^i by another complete orthonormal system of L2. It is well known (see, e.g., pQ]) that an 
operator F e £{^2) is Hilbert-Schmidt iff there exists a measurable function g : D x D ^ 'R such that 
F[u](-) — g{-,y) v{y) dy for v G L'^{D), and then, it holds that 

' dJ D 



(2.5) ||F||hs = / / g'ix,y)dxdy 



Let £hs(L2) be the set of Hilbert Schmidt operators of >C(L^) and <i> : [0,T] — > £33(^2)- Also, for 
a random variable X, let K[X] be its expected value, i.e., K[X] := J^_^XdP. Then, the Ito isometry 
property for stochastic integrals, which we will use often in the paper, reads 



(2.6) 
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For later use, we introduce the projection operator H : L {{0,T) x D) ^ L {{0,T) x D) defined by 
(2.7) n(g;.)|^^^ /" g{t,x)dtdx, VneM,, Vfiejf, 

for g e L^((0,T) x _D), which has the following property: 

4 



Lemma 2.1. For g G L^{{0,T) x D), it holds that 

(2.8) / / U{g;s,y)dWis,y)= f [ W{t,x) g{t,x) dtdx. 

Jo J D Jo J D 



Proof. To obtain (|2.8p we work, using ()2.7p and the properties of W, as follows: 

n n{g;s,y)dW{s,y)^^^ E E ( / 9 dtdx) ( T [ Xs^Js,y)dW{s,y) 
Jo Jo ^eA/'.^GJ," -^^"-^ -^o-^" 



-Atthy E E Cf 9it,x)Xs^Jt,x)R-'^dtdx 



Jo 



g{t, x) W{t, x) dtdx. 

□ 

We close this section, by stating some asymptotic bounds for series that will often appear in the rest 
of the paper and for a proof of them we refer the reader to Appendix [X] and Appendix [Bl 

Lemma 2.2. Let d G {1,2,3} and > 0. Then, there exists a constant C > that depends on and 
d, such that 

(2.9) HJ'''"'^ < C^-' Vee(0,2]. 

Lemma 2.3. Let d G {2,3} and S > 0. Then there exists a constant C > which is independent of d, 
.such that 

(2.10) J2 < c p.iSi) S"^, 
where pd{s) := 1 + J2i=i ■s'- 

2.2. Linear elliptic and parabolic operators. For given / G L'^{D) let Ve G H^{D) n Hq{D) be the 
solution of the boundary value problem 

(2.11) Ave = / in D, 

and Te : L^{D) — > H^{D) n Hq{D) be its solution operator, i.e. Tg/ :— Ve, which has the property 

(2.12) \\TEfL.,o< CE\\f\\„.-,,o: yf eH"'^-i"^^^'-'^D), ymeNo. 

Also, for / G L^{D) let Vg G H'^{D) be the solution of the following biharmonic boundary value problem 

(2.13) A^Ve = f in £>, 

(2.14) A'-Ve\^^^0, m-0,1, 

and Tg : L^{D) — > H''(_D) be the solution operator of (|2.13p . i.e. Tgf :— Vg, which satisfies 

(2.15) ||r,/|U.,z, < V/Gi?"'^'^^°'™-^>(i?), VmGNo. 
Due to the type of boundary conditions of (|2.13p , we conclude that 

(2.16) TEf = Tlf, yfeL\D), 
which, easily, yields 

(2.17) iTBVi,V2)o.D ^ {TEVi,TEV2)a.o y vi,V2 e L'^{D). 



Letting (iS(t)it;o)t6[o,T] be the standard semigroup notation for the solution w of (|1.3p . we can easily 
establish the following property (see, e.g., [23], [H]): for £ £ No, P, p £ and q e [0,p + 4£] there exists 
a constant C > such that: 

(2.18) j\T-taf\\dfS{T)wo\\l^dT<C\\wo\\l,^.,-2,-2 yh>ta>0, ^ Wo £ 11^+^'-^^-^ (D) . 

2.3. Discrete spaces and operators. For r £ {2,3,4}, we consider a finite element space Mh C 
Hq[D) n H^{D) consisting of functions which are piecewise polynomials over a partition of D in tri- 
angles or rectangulars with maximum mesh-length h. We assume that the space Mh has the following 
approximation property 

(2.19) inf \\v - xlU o < C h^-^ \\v\U, ^ yv £ H'^+\D)nH^{D), 

which covers several classes of finite element spaces, for example the tensor products of splines, 
the Argyris triangle elements, the Hsieh-Clough- Tocher triangle elements and the Bell triangle (cf. [7], 

A finite element approximation Vbm £ of the solution Vb of (j2.13p is defined by the requirement 

(2.20) BhVB.^ = Phf, 

and we denote by Tb.h ■ L'^{D) Mh the solution operator of (|2.20p . i.e. Tb,,,/ :— Vbm ~ Bf^^Phf for 
/ £ L^{D). It is easy to verify that ^ is selfadjoint, i.e., 

(2.21) {Ts.J,9)o.n = (/,T,,.5)o,z, V/,.g e L^{D). 

Also, using (|2?20ll . ([2?T3l) and ((2TT5| we conclude that 

||Ar«,J||o,„ <||AT^/|1„,„ 
(2-22) „ „ „ , 

<C||/ll-2,« yf£L^{D). 

Applying the standard theory of the finite element method (see, e.g., [7], [3]) and using (|2.15p . we get 

(2.23) ||A(T,/-T,,,J)||o,„< C/i'^-i 11/11.-3,1,, V/eiJ'°'^'^^'-3.o}(i?), 
while error estimates in the L^{D) norm are obtained in the proposition below. 
Proposition 2.1. Let r £ {2,3,4}. Then, it holds that: 

(h'Wfh^n, r = 4 

(2.24) \\Tsf-TB.J\\o.o<clh^\\f\\„^n, r = 3, y f £ H'^-^^^^-^-^^D). 

[h^\\f\\-,.n, r = 2, 

Proof. Let / £ ijmax{o,r-3}(^-) g^^^ g ^ _ y^^j ^Iso, wc define a bilinear form 7 : H^{D) x 
i?2(_D) ^ R by 7(wi,W2) := (Awi, Aw2)o,d for vi, V2 £ H^{D). Now, let w^, Ws £ 'H.'^{D) be defined by 
TflAe — Wa and T^e — Wb- Then, using Galerkin orthogonality, we have: 

l|Ve||^^ = - 7(u'A,e)o,i, 
'■^■^^^ <||Ae||o.i, inf Hwa-xIUd 

and 

l|e|lo,i, = -7(wA,e)o.D 

<||Ae||o,D inf - xI^d- 

Case i: Let r £ {2,3}. Then, using ([2^, ((2:23)) . ((2T9)) and ((2:22)) . we obtain 

||eL%<Cr-i||/||.-3,.,/i'-MhB||.+i,., 
<C/i2('-i)||/||^_3,„||e|j_3,, 

which, obviously, yields ()2.24p . 



(2.27) 



Case 2: Let r = 4. Then, combining, (EUl), ^^\^ and (HI]), we get 

|ie||2„ <C||Ae||o,„/i3||T«e||5.^ 
<C||Ae||„,„/i3||Ve||o,^. 

Also, we observe that (j2.25p and (|2.15p yield 

, , llVello.z, <||Ae|||„||r«Ae||L 

(2.28) 

<C||Ae||J^||e||?,„. 
Now, we combine (fS^ . ([2^ and (f2?23l) to have 

IjellL <C/i3||Ae|||.^ 

<C^/i'* ll/IlL: 

which obviously leads to (|2.24p for r 4. □ 

3. An Estimate for the Modeling Error 
Here, we derive an L^(Lp(L^)) bomid for the modeling error u — m, in terms of At and Ax. 

Theorem 3.1. Letu andu be defined, respectively, by (jl.2p and (|1.7p . Then, there exists a real constant 
C > 0, independent ofT, At and Ax, such that 



(3.1) max{E < C ^{pd{Ati))^ At^ + e'^ Ax—' 



Vee (0, ^1 



Proof. Using (|1.2p and (|1.7p . we conclude that 

(3.2) u{t,x)~u{t,x)^ f f [X^n,t){s)G{t-s;x,y)-G{t,x;s,y)]dW{s,y) V {t,x) e [0,T] xD, 

Jo J D 

where G : (0,T) x D ^ i^((0,T) x Z?) given by 



G(t,x;-) 



At (Ax) 



X(o,t){s')G{t-s'-x,y')ds'dy', ^neK. V^* e J^. 



Let e := {E [\\ul - Utllo ^] }' and t e (0,T]. Using ((^ and Ito isomctry we obtain 



0W = 



At (Ax)'' 



E E 



'^(o,t)(s)G(t- s;a;,y) 



A'(o,t)(s')G(t-s';x,2/') ds'dy 



dsdy > 



Now, we introduce the splitting 

(3.3) 

where 



e(i) < e^(t) + e,(t). 



e.4(i) := 



At (Ax)'' ' 



E E 



'^(o,t)(s) G{t-s;x,y) 
— G{t — s;x,y') ds'dy' 



dsdy > 



and 



Xio.t){s) G{t - s;x,y') 



— A'(o,f)(s') G{t — s'] X, y') ds'dy' dsdy > dx 



Estimation of 0a(O- Using (jl.Sp and the (•, orthogonality of {ea)a£nd, we have 



-^(o.t) (s) G{t- s;x,y)-G{t- s;x,y') dy' 



{Ax 



r)2d E E 



= (A^)^E IE / A'(o,,)(s)e-2^=(*-^'dsU ^ / (/ (e„(y)-£„(y'))rfy')'d2/ 



from which, using the Cauchy-Schwarz inequahty, follows that 

(3.4) elit) < E ( ^-"^"-^^ ^■^) f (A^ E / 



«(y)-ea(y')| c^y'c^y 



Observing that e"^^"^*"''^ ds < i for a e N'^, and that 

sup \ea{y) -£a{y')\ <25+i min|l,f d3 Aa;|a|„d| 
v,v'&D^ J 

<2l+i-7^7rfi Aa;T|a|^,, V7e[0,l], Va G V/i G J-/, 

(03) yields 

(3.5) ei(t) < 2^+1-2-^ TT^^-^ (Aa;)^^ E 

The series in p.Sp converges when 2(2 — 7) > d or equivalently 7 < Thus, combining p.Sp and 

(|2.9p . we, finally, conclude that 

(3.6) eA(t) < Ce-^ Ax^-' Vee(0,^]. 



Estimation of QB{t): For t E (0,T], let iV(i) min { £ G N : 1 < ^ < and t < } and 

f r„, if n < Nit) 

T„(i) :=T„n(0,t) = <^ ^ , n=l,...,N{t). 

[ (i^(*,-i,t), if n = iV(t) 

Now, we use (jl.Sp and the (•, •)o o— orthogonality of (eq.) follows 



(Ax)" 
(At(Aa:)<i)2 



E E 



(At (Ax)")^ 



E 



X^oM^) G{t - s; y') - A'(o^t)(s') G(t -s';x, y') 

N(t) 



ds'dy' 



ds 



E(/ ^M')dy'y E/(/ (A'(o,)(.)e-^=(*-^)-A'(o,)(s')e-^"(*-^')) 



which yields that 

/ N{t) 

(3.7) e|(i)< 2<^ IMP E*nW 1 , 



where 



Kit)--= I ( I (A'(o,t)(s)e-^^(*-«)-A'(o,t)(s')e-^^(*-'*')) ds^)' ds. 
Let a e N'^ and n G {1, . . . , N{t) - 1}. Then, we have 
( If Xle-^l^'-^Urds'Yds 



p , p /'max{s',s} 2 

< U Xle-'^"o.(t-r) drds') ds 



< 2 



L U^It ^^c.e-^"^'~^^dTds'yds + 2j^ (^j^f Xle-^'^^'-^Urds'y ds 

from which, using the Cauchy-Schwarz inequaUty and integrating by parts, we obtain 

nit) < 4 (At)2 J (e-^'(*-«) - e-^'(*-*"-i) )' ds 

< 4(Ai)2(l-e-^«^*)' / e-2^"(*-«)ds 

< 2 (At)2 (1 - e-^>*)' ^'-'"'-'^"''-'"-^ - 



A2 

Thus, by summing with respect to n, we obtain 

w(t)-i 

(3-8) ^ E *nW< 2 XT 

Considering, now, the case n = N{t), we have 

(3.9) n,.,it) = n{t) + mt) 

with 

ft 



^1{t)-= f if rA2e-^«(*-)dTds'+ r'""e-^'(*-«)ds') ds 



-N{t)-1 \ -N{t)--L 



n{t):= I ( / e-^i(*-')ds' ) ds. 

Then, we have 



*iV(t) ^ 



1 — e~^i (*-%(*)-! ) 



and 



N(t)-i L ^/v(t)-l 



ds 



< 2 



< 2 



JV(f)-l 

t /"maxlsjs'} 



ds + 



A?, 



1 _ g-2Ao(*-t]v(t)-i ) 



W{t)-1 W(t)-1 



N(t)-1 
t 



< 8(At)2 

< 8(At)2 

'■?f"(t)-i 

which, along with (13. 9|) . gives 



+ ( 1 - e 



(l_e-2^=^*) 



+ ( 1 _ e 



< 5 ( 1 - e-^= ) + t ( 1 - • 
Since the mean value theorem yields: 1 — e^'^^'^* < A^ At, the above inequality takes the form 
(3.10) 

Combining ^1}, and ([XTU]) we obtain 

(3.11) elit) < 8 ^ 

Now, combine ((3lT|) and ((TTO)) to arrive at 

(3.12) e^(t) < c(pd(Ai*))5 Ar 

The error bound p.ip follows by observing that 9(0) = and combining the bounds p.3p . (|3.6p and 

4. Time-Discrete Approximations 

The Backward Eulcr time-discrete approximations to the solution u{Tm,-) of the problem (jl.6p are 
defined as follows: first, set 

(4.1) {7° := 0, 

and then, for m 1, . . . , M, find [/" e H'*(i:)) such that 



1 , . 1 . 4-d 



(4.2) 



2 /- j m 



W ds a.s. 



To develop an error estimate in a discrete in time L'^{Lj,{L'l)) norm for the above time-discrete ap- 
proximations, we need an error estimate for the Backward Euler time-discrete approximations, (W"^)^^q, 
of the solution w to the deterministic problem (|1.3|) . given below: First, set 

(4.3) W° := wq. 
Then, for m = 1, . . . , M, find VF" G ii'^{D) such that 

(4.4) W"-M/™-i-f fc„,A2w^™ = 0. 
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Proposition 4.1. Let (VF™)^^^q be the Backward Euler time- discrete approximations oj the solution w of 
the problem (jl.Sp defined in (|4.3p - (|4.4p . If Wq G H^(I?), t/ieTi, there exists a constant C > 0, independent 
of T , At, Ax, M and {km)m^i, such that 

(4.5) ( Efcm||W^"-i«(T™,-)llL ) < C^(^»aj'lko||H«- V0e[O,l]. 

\ m=l / 

Proof. The proof is omitted since it is moving along the hncs of the proof of the one dimensional case 
which is exposed in Proposition 5.1 of p[8]. □ 

Next theorem proves a discrete in time L^(Lp(L^)) convergence estimate for the Backward Euler time 
discrete approximations of u, over a uniform partition of [0,T]. 

Theorem 4.2. Let u be the solution of (|1.6|) and {C^"'}m=o Backward Euler time-discrete ap- 

proximations specified in (|4.ip - (|4.2p . // km — Ar for to = 1, ... , M , then there exists constant C > 0, 
independent of T , At, Ax and Ar, suc/i that 

(4.6) max |E[||C/™-u(r„,.)||Lir < ^(AT,e) At^-% Vee(0,^], 

l<m<i\/ I L ' J J " 

w/iere w(At, e) := [e"5 + (At)^ (^^(Ati))^]. 

Proof. Let / : ^^(i:)) L^{D) be the identity operator, A : L'^{D) H^(D) be the inverse elliptic 
operator A := (/ + ArA^)^-^ which has Green function Gf^{x,y) = X^agN'' ^i+AtA''^^ ' ^-^^ ^/(^) ~ 
J^GAix,y)f{y)dy for x G D and / € L'^{D). Obviously, Gf,{x,y) — Gf,{y,x) for x,y £ D, and G G 
i^(_D X _D). Also, for to G N, we denote by Ga,™ the Green function of A™. Thus, from (|4.2p . using 
an induction argument, we conclude that C/™ = X^jLi /a A™^-'+^W^(t, •) dr for to = 1, ... , Af , which is 
written, equivalently, as follows: 



(4.7) U"'{x)= / JCra{T-x,y)W{T,y)dydT Vx G A to=1,...,M, 

Jo J D 

^here JCm{T\x,y) ■.^YJjLi^/^j{T)G,,^m-]+i{x,y) VTG[0,r], ^x^yeD. 

Let TO G {1,.. .,Af} and := E[||C/" - w(t„, .)||2^] . First, we use dUIl), (dUl), ([211) and ([231), to 
obtain 

= 



'^(o.r„.)(T) [^m(T;x,y) - G(t„ - t;x,?;)] VK(t, y) dydr^ dx 
/ S X! X! / '^(o,r„)(T") [^Tn(T;x,2/) - G(t™ - r;x,y)] drdy Wx 



At (Ax)'' 

Then, we apply the Cauchy-Schwarz inequality and p.5p to arrive at 



f ™ < y" (^J J [i^m{T;x,y) - G{T,n - T;x,y)]^ dydx^ dr 

51/ ( / / [GA,m-e+iix,y) - G{Trn - T;x,y)]'^ dydxj dr. 

J Aj \J dJd J 



Ihs ' 



< 



<^ / ||A"-^+i-5(r„,- 
Now, we introduce the splitting 

(4.8) < B'^ + 6^, 
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where 



m „ 

Sr:=2^/ ||A 

m „ 



m-e+1 



SiTm~Ti^l)\\lsdT, 



T,n - - 5(r,„ - r)||^ dr. 



Estimation of S™: By the definition of ttie Hiibert-Sclimidt norm, we have 

m / 



=1 Xaer 



< 2 Ar||A™-^+i£,-5(T™-r,_i)efe||2„ 

< 2 5] [f] Ar||A^efe-5(r,)efe||2 J. 
Let 8 £ [0, ^-^)- Using the deterministic error estimate (|4.5p . we obtain 



(4.9) 



The convergence of the series is ensured because 4(1 — 29) > d. 

Estimation of B^: Using, again, the definition of the Hilbert-Schmidt norm we have 



(4.10) 



Observing that S{t)ea = e~^°* Sa for t > 0, (liTTO]) yields 



B^=2J2 
= 2E 



=1 



E 

m 

E 



£Q,(a;) da; I dr 



<2 E (1-^ 
< E^ 



= 1 " 

-a! At \2 



1 _ g-^aC'^— ^f-i) 
e-2A^(r„-r)^^ 



dr 



A?, 



from which, applying (|2.10p . we obtain 

(4.11) B^ < C PdiAri) Ar^. 

Thus, we obtain the estimate (|4.6p as a conclusion of (|4.8p . (|4.9p . (|4.1ip and 
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5. Space-Discrete Approximations 

Let r e {2, 3, 4}. The space-discrete approximation of the solution u of (|1.6p is a stochastic function 
Uh ■■ [0, T] Mh such that 

dtUh + BhUh = PhW on (0, T] , 
5.1) a.s.. 

u,m = 

To develop an L"^ {Lp{L^)) convergence estimate for the space-discrete approximation Uh, we will derive 
first an L^^L"^) error estimate for the corresponding space-discrete approximation Wh of the solution w 
of lOl) (cf. [llj and [2]), which is a function wh ■ [0,T] ^ M,, such that 

dtWh + BnWh^Q on (0,T], 

Wh(0) = PhWQ. 

Since Wh can be considered as the value of a linear operator of the initial condition wq, we will write it 
as Wh{t, ■) — [Sh{t)wQ\{-) for t £ [0,T]. Thus, by Duhamel's principle (cf. [13]), we have 

(5.3) Uh{t,x) = f [Sh{t- s)W{s,-)]{x)ds a.s.. 



Proposition 5.1. Let r G {2,3,4}, w be the solution of (jl.Sp and Wh be its space- discrete approximation 
given in (j5.2p . If Wq £ H^(-D), then, there exists a constant C > 0, independent of T and h, such that 

(5.4) (^J^ \\w-Wh\\l„dt^ < C h""^-''^ \\wo\U^.e, yOe[0,l], 

where 

{26 ifr = 2 r36'-2ifr = 2 

46 if r = 3 and ^{r, 6) := I 49 - 2 if r 3 . 
5 6 ifr = 4 [56'-2ifr = 4 

Proof. Let e := w — Wh and p :— [Tb^u — Tb)A'^w. We will derive (|5.4p by interpolation, after showing 
that it holds for 6* = 1 and 6* = 0. 

Observing that T^jiCt + e = p on [0, T], and then taking the (•, •)o,d inner product with e, we easily 
arrive at 



(5.6) / \\e\\l^dt< / WpWl^dt. 

Jo Jo 

For r = 2, using (g^, fTM^ . ([^ and (^JE\\ . we have 

(5.7) (f l|e|lL^t)^<C/.^(/^HI^3d. 

< C/l2||«;o||^i. 

Also, for r = 3,4, combining dH]), (PTM]) . (jO)) and ((TT5)) we get 

(5.8) (f l|e||^,..^t)^<C/.-(/^|lH!^..rf.' 

< c;i'-+i|kolU^-i. 

Thus, relations (|5J)l and (|5?8)l yield ((O)) for 6* = 1. 

Since Tsiyf + w = on [0, T], we obtain (T^wt, w)o,d + ||w||n „ = on [0, T], which, along with (|2.17p . 
yields ^llT^wllj^ o + = on [0, T]. Then, integrating over [0,T] and using ((2T2)) . we get 

(5.9) (^^"||Hlo%'^^)' < C|koll-.,o. 
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(5.10) 



Since Te^dtWh + Wh on [0, T], we obtain (TB,hdtWh,Wh)o,D + \\wh\\o^o =0 i^^T], which, along with 
(I2:2TD . yields ^(T^.^^w/., ^y/^)o,D + \\wh\\lr, = on [0,T]. Then, integrating over [0, T] and using (E^H), 
we have 

< ||ArB,,,Wo||o,D 

< C \\wq\\_2.D- 

Hence, from (|53)) . ((5?T0)) and (|2^ . we obtain (/q^ ||e||^ „ di) ' < C||wo||h-2, which yields with 
6' = 0. ' □ 

Next lemma shows that a discrete analogue of ()1.4p holds. 

Lemma 5.1. Let r e {2,3,4} and Wh be the space- discrete approximation of the solution w of (jl.3p 
defined in ()5.2p . Then, there exists a map Gh ■ [0, T] — > C(I? x 13) smc/i that 

(5.11) Whit;x) ^ [ Gh{t;x,y)wo{y)dy Vte[0,r], Vx e A 

and G'/i(t; a;, y) — Gh{t; y, a;) for x,y £ D and t £ [0, T] . 

Proof. Let dim(M/i) — Uh and 7/; : Mh x Af^ ^ R be an inner product on given by 'yhiXAiXs) '■— 
{Axa, AXb)o,d for Xa, Xb G Mh. We can construct a basis (Xi)j=i of Af/i which is iy^(_D)— orthonormal, 
i.e., (xijXi)o,D = (^ij for i, j = 1, . . . ,nh, and 7/1— orthogonal, i.e., there are {\h,tyi!Li C (0, +00) such that 
IhiXi: Xj) — ^h.i Sij for i,j — 1, . . . ,nii (see Section 8.7 in [9J). Thus, there exists a map uj : [0, T] R"'' 
such that Wh{t]x) — J^^^^i'^ji^) Xj{^)- Since Wh{0) = PhWq, it follows that Wj(0) = (wo,Xj)o,D for 
j = l,...,n/i. Now, dSSl) yields that J^w(t) = Buj{t) for t e [0,r], where B £ M"'^x"'> with B^j := 
-7?i(Xi: Xi) = ->^h,i Sij for j, j = 1, . . . , n/i. Hence, it follows that Lu^{t) = e'^'^'" * (wq, X/!)o.d for t £ [0, T] 
and £ = 1, . . . ,n/,, which yields (|5.1ip with Gh{t]X,y) = Yljti (^'^"-^ *Xjix)xj{y)- □ 

We are ready to derive a convergence estimate, in an L'^{L'^{L1)) norm, for the space-discrete ap- 
proximation Uh to the solution u of the regularized problem. 

Theorem 5.2. Let r £ {2,3,4}, u be the solution of p.6p and Uh be its space- discrete approximation 
defined in (|5.ip . Then, there exist a constant C > 0, independent of T , At, Ax and h, such that 



(5.12) max{E[||u,,-^2||„%]}^ < Ce-U''(''''^), Ve G (0, K^, d)], 

where 

f ^ if r 2 



^ if r==3,4 



(5.13) v{r,d): 
Proof Let e:^Uh~u and t £ (0, T]. Then, (|5T11) and dTT]) yield 

e~(t,a;)= / / ^Gh{t ~ s]x,y) ~ G{t — s;x,y)^W{s,y)dsdy VxeD, a.s.. 
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Thus, using the ltd isometry property of the stochastic integral and the Cauchy-Schwarz inequahty, we 
have 



E[|le(t,.)IL%] =E 



D \J0 Jd 



At (Ax)" 



< 



lO \JDJn 

which, along with (|2.5p . yields 
(5.14) 



Xio.t){s) [Gh{t ~ s;x,y) ~ G{t - s;x,y)\ W{s,y) dsdyj dx 

Xio,t){^') [Gh{t - s'; X, y') - G{t - s'; x, y')] ds'dy'^ dx 
[Gh{t - s; x, y) - G{t - s;x,y)Y dydx ) ds 



E[l|e(i,-)IIL] < / ||5(,s)-5,(s)||;^ds. 



Since e(0, •) = 0, we use (I5.14p . the definition of the Hilbert-Schmidt norm and (|5.4p . to obtain 
maxE[|le||y < / ( ^ || 5(s)e„ - 5/,(s)e„|L^^ ) 
- E ( / Pi^)^k-Shis)ek\\l^ds 

(5.15) 



2«(r,e) 



The series in the right hand side of (|5.15p converges if and only if — 2^(r, 6) > d. Thus, in view of (|2.9p . 
we arrive at ((5TT2| and ((57T3)l . □ 



6. Convergence of the Fully-Discrete Approximations 
6.1. Consistency estimates. First, we derive some Holder-type bounds for u. 



Lemma 6.1. Let u he the solution of (|1.6p . Then, there exist a real positive constant C, which is 
independent of T , At and /S.x, such that 



(6.1) 

and 
(6.2) 



E 



2(t6, •) - u(t, •)] dr 



< C {Pd{{n - Ta)-^))- \n-Ta\'+- 



for Ta, Tb e [0, T] with Ta < Tb. 

Proof. We will omit the proof of (|6.2p because it is similar to that of (|6.ip which follows. 

Let m € {1, ... , M}, n e (0, T] and Tq e [0, T] with Ta < n, and ^(•) := JJ^ [uin, •) - w(r, •)] dr. 
First we assume that Tq > 0. Then, we use (|1.7p . (|2.8p . the Ito- isometry property of the stochastic 
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integral, (|1.5p and the L^(i))— orthogonality of i£a)ai£N''i to obtain 



E E[ 



[^('a,Tt){s')G{n ~ s';x,y') 



Xio,r){'')GiT-s';x,y')] drds'dy' 



dx 



E E E 



['^(o,r,)(s )e 



which, along with the use of the Cauchy-Schwarz inequahty, yields 



ie[IImIIL] < E 



'^(o,r,)(s )e 



'^(o,r)(s )e 



(6.3) 



< in - Ta) J2 



ds'dr 



-^^l(rb-^')ds'dT 



< 



„-2 A;i (t^-t 



Finally, we combine (j6.3|) and ()2.10|) to arrive at ()6.1|) . The case Ta — follows by moving along the lines 
of the proof above using that u(0, x) — 0. □ 

Next, we show a consistency result for the Backward Euler time-discrete approximations of u, which 
is based on the result of Lemma 16.11 

Proposition 6.1. Let u be the solution of (|1.6p and {^m)m=i stochastic functions defined by 



(6.4) 



u{Tmr) -u{Tm-ir) + hn^ u{t„i,-) = I W dr + a.s., m = l 



,M. 



Then it holds that 

(6.5) {E[\\T^d^\\l^])^ < C (pdikk))-' {k„,)^+^, m = l,...,M. 

Proof. Let m G {1, . . . , M}. Integrating the equation in (|1.6p over and subtracting it from (|6.4p . we 
conclude that Tsdm{ ) = [^('''mO ~ •) ] dr a.s.. Thus, to get the bound (|6.5p . we apply the result 
(|6.ip on the latter equality. □ 



6.2. Discrete in time (Lp(L^)) error estimate. We first obtain a discrete in time (L^(L^)) error 
estimate for the Backward Euler fully-discrete approximations of u, by connecting it to the error estimate 
of Theorem 15. 21 for the space-discrete approximation of u. 

Theorem 6.2. Let r e {2,3,4}, u be the solution of ([T^ and 0^)'^^=^ C Mu be the Backward 
Euler fully-discrete approximations of u defined in (jl.lOp - p.lip . Then there exists a constant C > 0, 
independent of T , r, Ai, Ax, h, M and {km)m=i! such that: 



(6.6) 



X^A:„E[||&r-^(T™,-)llLl < CVf\e~^- h'^^^^"'^"' +w{k^^^) (fc^)" 



for e & {0,i^{r,d)], where Lj{k,„^^) :— (pd((^max) * )) ^ andv{r,d) is defined in (|5.13p . 
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Proof. Let Uh be the space-discrete approximation of u defined in (jS.ip . e — u — Uh, z'^ '■— UJ^^ ~Uh{Tm) G 
S";; for m = 0, . . . , M, and Vh := {Em=i [H^rllo.i,] First, we observe that 



(6.7) 



^fc„,E[||&™-u(r„„.)|l„% 



< Vh + VT ma^{E[\\e\\ljy 



Integrating ()5.1|) over Am and subtracting the obtained relation from we arrive at 

(6.8) r«.,(zr-4""') + fcm^r =P/>.rn a.s., m=l,...,A/, 

where Ph,m ■— [uh{T, ■) — uii{t„i, •)] dr. Take the (•, •)(,o~ inner product of both sides of (|6.8p with 
z™, sum with respect to m from 1 up to M, to obtain 



(6.9) 



mil 2 



h.m.-) )q,d ^-S. 



Since zO = 0, we conclude that E;^^i(Ar«,,z™- AT^,,z™-i, Ar«.,z,7)o,,, > i || AT, 
taking expected values in (|6.9p and using the Cauchy-Schwarz inequality we get 



a.s.. Thus, 



(Vh)^ <E 



m— 1 



(6.10) 



< E 



E 



_i J D J A 



M „ 

< ^ / E[P„(T,.)-Uh(rm,-)IIL] rf^- 



Using (|6.10p and (|6.2p . we conclude that 
(6.11) 



Vh<\j2 I E[lle(T,.)-e(rm,-)llU dr 

11) lm=l-^'^™ 



E^ E\\\u{t,-)-u{t„„-)\\1^] dr 



<CVT 



max I 

[0,T] 



+ (pd((fc._)'))^ (fc._) — 



Thus, (HH) follows from ([131), ([OT|) and ((51^ . 



□ 



6.3. Discrete in time L^(_Lp(_L^)) error estimate. To get a discrete in time i^(Lp(L^)) error 
estimate for the Backward Euler fully-discrete approximations of u, we compare them to the Backward 
Euler time-discrete approximations of u defined in (|4.ip - (|4.2p . 

For that we derive first a discrete in time L^{L'^) error estimate between the Backward Euler time- 
discrete and the Backward Euler fully discrete approximations of the solution w of (|1.3p given below: 
First set 

(6.12) := PhWo. 
Then, for m = 1, . . . , M, find W^^" G Mh such that 

(6.13) - W;^'' + BhWl^ = 0. 

Proposition 6.3. Let r e {2,3,4}, w he the solution of the problem (|1.3p . (M^™)m=o Backward 
Euler time-discrete approximations of w defined in (j4.3p - (l4.4p . and (W^™)m=o Backward Euler 

fully-discrete approximations of w specified in ()6.12p - ()6.13p . //wq G H'^(I?), t/ien, t/iere exists a constant 
C > 0, independent of T , h, M and (fc„i)m=i7 si^c/i that 



(6.14) 



5] fcm liw^" - wnlo < /^"^^^''^ lkolU,>,«, V0 e [0, 1] 
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where v{r,9) and ^{r,6) are defined in (j5.5p . 

Proof. Let E"' := W"' - WJJ" for m = 0, . . . , M. We will get by interpolation, showing it for 6* = 

and 9^1. 

We use g31) and (O^ . to obtain: T^^uiE"" - E""-^) + E"' = fc™ (T^ - T^.^A^iy™ for m = 
1, . . . , M. Since Tg j^E'^ — 0, proceeding as in the proof of Theorem 16. 2[ it follows that 

M M 

(6.15) J2 ^™ ll^^llo... < E ^™ IK^- - T,m)A'W"^\\1^ . 

m— 1 m— 1 

Let r = 3. Then, by (P?^ and (pT5)) . we obtain 

M M 

(6.16) E k^WE^lln < Ch' ^"^ ll^'^"!!'.- 

m— 1 rn — 1 

Taking the (., .)o,r)— inner product of (|4.4p with A^M^™, and then integrating by parts and summing with 
respect to m from 1 up to A/, it follows that 



(6.17) E(AH^™- AH^™-\AH^™)o,o + E fc™ || A^M^™!!^ ^ = 0. 

■m— 1 m— 1 

Since E™=i (AW" " AW^"-i, AT^")^^ > i ( |1AW^-||2^ - \\AW°\\ln), yields 

(6.18) EA;™||A2M^"||2^< i||u;o||?,.. 

m— 1 

Combining, now, (|6.16p . (|6.18p and (|2.3p . we obtain 



(6.19) 



E ll^™llo,x. < C/i^llT^ollH- 



Let r = 2. Then, by jlJll, (El) and (|6T5l) . we obtain 



(6.20) 



m— 1 

A-/ 

-^fc™ (AW^™,A2w^'")o,^ 



m—1 



Taking the (■, ■)o o— inner product of (|4.4p with AVl^"^, integrating by parts and summing with respect 
to m from 1 up to Af , it follows that 

M M 

(6.21) (VM^'" - VM^™-\ VM^")^ z, - ^™ (A^W^", AWK'")o,,, = 0. 

m—1 m—1 

Since E™=i(VVF™ - yW"'-\VW"')o,o > H II^W^"llo,i, " l|VW^°||2^], (EUD yields 

A/ 

(6.22) - Y (A^W^", AW^™)„,„ < i \\w4l^. 

m—1 

Combining (P?^ and (US]) we get 



(6.23) 



Y^raWE^Wlo] <Ce\\w,\W. 
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Let r = 4. Then, observing that A'^W"^ G ii^{D) and using the relations ((^ and (pTTC)) . we 

obtain 



m— 1 



(6.24) 



m— 1 

M 

<Ch^° km \\A^W'^\\l_, 

M 



m— 1 
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m— 1 



After, Applying the operator A on (|4.4p . take the (•, •)o.£,— inner product of the obtained relation with 
A'^W"\ integrate by parts and sum with respect to m from 1 up to AI, to get 



(6.25) 



Also, we have 



(6.26) 

Thus, fOS)) and fOel yield 
(6.27) 



J2 (AT^" - AM^™-\ A2m^")o,b > J2 ( II AW^"llfti - I|AT^"|1hi IIAiy 



7n — l\ 



>i (||AW^-||2, -||AT4^°||^/ 



- J2 km iA^W"^,A'W"^),,^ < \ ||u;o|l|3 



Combining (pTM)) and we get 



(6.28) 



E 11^' 



m||2 
O.D 



< C}t'\\'W^\\^^. 



^ m— 1 



Thus, the relations (jO^ . and (j^:^ yield (jO^ for = 1. 

Since Tb,/,(W^^ - W^'^) + fc,„ VK^ = for m ^ 1, . . . , M, we obtain 

(AT,,,t^,r - AT^,,W^-\AT^,^M]^\^r. + fc™ ||W^rilo,z,, m = 1, . . . ,M, 
which, along with and ([^ . yields 



(6.29) 



E^™ii^"iio%<liiAr«.'^^oii^,. 



m— 1 



< C ||wo||h-2. 



Now, using gll) and (l2l7l) . we obtain (T^M^" - TsVK™-\ TeVK™)o,b + fc™ ||W^™||o,o = for m = 
1,...,M, which yields - WT^W^-^H,, + 2 fc™ ||T^™||2 ^ < for m = l,...,Af. Then, 

summing with respect to m from 1 up to Af , and using (|2.12|) and (|2.4[) we obtain 

A/ 

E^-iiw^"iio,.< ir^"'°iio% 

< C||u;o||^-.. 



(6.30) 



fe=i 



Finally, combine ([QQ]) with (lOO]) to get ( Em=i ^™ ll^™llo d ) ' < C'lkolU-i, which is equivalent to 
(|6T4l) for 61 = 0. □ 
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The following lemma ensures the existence of a continuous Green function for the solution operator of 
a discrete elliptic problem. 

Lemma 6.2. Let r e {2,3,4}, e > 0, / G L'^{D) and iph e Mh such that 
(6.31) eBni^h+^Jh = Phf. 



Then there exists a function Gh^e G C{D x D) such that 

(6.32) M^)^ I GH4x,y)f{y)dy VxeD 

J D 

and Gh.e{x, y) = Gh,eiy, x) for x,y e D. 

Proof. Keeping the notation and the constructions of the proof of Lemma 15.11 we conclude that there 
are (Aij)"=i C M such that iph = YTjtii^jX]- Thus, (|6.3ip is equivalent to = i+lxu ^ (f^Xi)o.D for 

i = l,...,nh. Finally, we obtain (I02l) with GhA^, y) = Y^-^i ^i+lll^f ■ ^ 

We are ready to compare, in the discrete in time LJ"(L^(L^)) norm, the time-discrete with the fully- 
discrete Backward Eulcr approximations of u. 

Proposition 6.4. Let r E {2,3,4}, u be the solution of the problem (|1.6p . (C/™)m=o Backward 
Euler fully- discrete approximations of u specified in (|1.10p - (jl.ll[) . and (?7'")m=o Backward Euler 

time-discrete approximations of u specified in (|4.H) - (|4.2p . If the partition {Tm)m=o uniform, i.e. k„i = 
Ar for m = 1, . . . , M , then, there exists a constant C > 0, independent of Ax, At, h, M and Ar, such 
that 

(6.33) max (En|[/;;"-C/"f j|' < Ce-^ \/ 1 & {0,v{r,d)] 

l<m<M L L 110, r>j J 

where h'{r,d) has been defined in (j5.13p . 

Proof Let / : L^{D) L^{D) be the identity operator and Ah : L'^{D) -> SJ^ be the inverse discrete 
elliptic operator given by ■= {I + At Bh)~^Ph and having a Green function Gh^&^ (cf. Lemma 
Also, for £ S N, we denote by Gii,&T,e the Green function of A^. Using, now, an induction argument, from 
(|l.lip we conclude that C/S™ = X^ILi /a A!""-'"''^ VF(r, •) dr, m = 1, . . . , M, which is written, equivalently, 
as follows: 

(6.34) Uj:'ix)= f "1 Vh,miT;x,y)W{T,y)dydT ^xeD, m=l,...,M, 
where 

m 

2?/i,m(T;a;,?;) := ^A:'A^.(T)G'/i,A^,m-j+i(a;,?;) VTe[0,r], yx,yeD. 
i=i 

Using (j4.7p . (|6.34p . the Ito-isometry property of the stochastic integral, (|2.5p and the Cauchy-Schwarz 
inequality, we get 



[£™(t; X, y) - Vh,m{T; x, y)] ^ dydx^ 



W'"" -U'nlo\< ^ I I \^m{T;x,y)-Vh.^(T-x,y)V dydx]dT 

m „ 

<^/ ||A™-^-+i-Ar^-+i|lLdr, m = l,...,M, 
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where A is the inverse eUiptic operator defined in the proof of Theorem 14.21 Now, we use the definition 
of the Hilbert-Schmidt norm and the deterministic error estimate (j6.14l) . to have 



E 



\U'"-U, 



ni 1 1 2 
h Wo.D 



EiiA 



m-j+l 



A' 



m~j + l 



aeti'' 



< C/i^^M) ^ ||e„||2-,^,,,, m = l,...,M, V0e[O,l] 



Thus, we arrive at 



max E 

l<m,<M 



a 



2«(r,e) 



from which, requiring ~2^{r,9) > d, (j6.33p . easily, follows (cf. Theorem 



V^e [0,1], 



□ 



The available error estimates allow us to conclude a discrete in time i^(Lp(i^)) convergence of the 
Backward Euler fully-discrete approximations of m, over a imiform partition of [0,T]. 

Theorem 6.5. Let r S {2,3,4}, v{r^d) he defined by (j5.13p . u be the solution of problem (|1.6p . and 
(^™)m=o Backward Euler fully- discrete approximations of u constructed by p.lOp - dl.lip . If the 

partition (Tm)m=o uniform, i.e., km = Ar for m = l,...,M, then, there exists a constant C > 0, 
independent of T , h, Ar, Ai and Ax, such that 



max E 

0<m<M 



[v.\^\U^^~u{t„,,-)\\1,]Y 



< c 



w(At, ei) At V-^i + ^ h"'^''''^^-'' 



/or eiG(0,^] anrfea e (0,z/(r,d)] w/iere w(Ar, ei) := ^ + (Ar)^i (prf(Ar3)) 
Proof. The estimate is a simple consequence of the error bounds ()6.33p and (|4.6f 



□ 
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Appendix A. 
Proof of ([^ . Let e e (0,2]. First, we observe that 



00 /--l-oo 



dx 



< (2+1 ^ 1 



which easily yields (|2.9p for d = 1. For d = 2, we have 

00 „ 

j: iai;.(^--) < 2 j: — + / d. 

<2V;^+/ I r-^'+'^'Urde 
For c? = 3, using (|2.9p for 0? = 2, we proceed similarly as follows 



QGN3 /3eN2 "'(l,+-oo)3 

<3 V + r [" [' smi9)r-(^+'*'Urded(j) 

a^KM9 Jl Jo Jo 



□ 
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Appendix B. 

Proof of sum . First, we recall from [H] that J2T=i ^^1^"^^^' ^ C{l + Si)di. For d = 2, using the 
latter inequality, we have 

— xi — - ^ ^H^+k'^y^ — ^ / — ^^^w — 

agN2 n=l "'(l,+oo)2 



oo 



n=l -^0 



/' + OC 



Jo 

I r+oo 

<C(l + 53),5i + |l / e-"'dz, 
Jq 

which yields X^aeN^ " — C (1 + 6^ + 5^) . Finally, when o? = 3, using (|2.10p for d = 2, we obtain 



aeN3 /3eN2 " -'(i.+oo)^ 



/O "'0 "'0 

Jo 

<C(l + 53+j3),5H|^ / e-"'dz, 

JO 

which yields J2c.em aV° ' - C (1 + 6i + + Si) Si . □ 
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